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Abstract 

We report the experimental observation of the recently predicted peaks on 

the I- V curve of the superconducting single-electron transistor at relatively 

high temperatures. The peaks are due to the matching of singularities in the 

quasiparticle density of states in two electrodes of a tunnel junction. The 

energy shift due to Coulomb blockade provides the matching at finite voltage. 
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SinEfle-electron effectJ in superconducting structures have several additional featuresB'i 
in comparison with that in normal metals or semiconductors. The main differences are due 
to the specific role of the parity of the electron number on a small island,i'@ the effects of the 
Josephson coupling,0'iH and the specific shape of the quasiparticle density of states (QDS). 
The last topic received relatively small attention so far in both theoretical and experimental 
single-electronics, although QDS leads to various interesting effects. Besides the well-known 
shift of the Coulomb blockade threshold by 4A/e in SSS single-electron transistor (SET) 
(by 2A/e in NSN or SNS cases), let us mention the direct reproduction of the QDS on the 
/- V curve of the SET with discrete energy spectrum of the central electrode! and in the 
case of odd-parity current,0 singularity-matching (SM) peaks at finite temperatures,^ the 
limitation of the differential conductance by quantum resistance0'0 and the conductance 
jump at V" = 4A/e due to cotunnelingS. 

In this paper we discuss the theory of SM peaks in more detail and report their experi- 
mental observation in the SSS SET. Somewhat similar experimental results will be presented 
soon by another group.El 

The origin of SM peaks can be easily understood from the "orthodox" theory of the 
SET.0'0 Let us neglect the effects due to Josephson coupling and consider only the quasi- 
particle tunneling. The dc current / through the SET consisting of two tunnel junctions in 
series can be calculated from the equations00 

/ = E[r^W-rrH]^H, (i) 

n 

a{n) [r+(n) + T^{n)\ = a{n + 1) [T^in + 1) + r+(n + 1)] , (2) 

where Tf(n) are the rates of tunneling through ith. junction (i = 1,2) in the positive (-I-) 
or negative (— ) direction when n extra electrons are present on the central electrode of the 
SET, and a{n) is the probability of the charge state n. In the case of SS junction the 
tunneling rates Ff (n) = Ti(Wj^{n)) are given by equationJl 

1 f + OO 

r^W = ^/ p{e)f{e)p{e + W)[l-f{e + W)]de, (3) 



where p{e) is the normahzed QDS, A(T) is the superconducting energy gap, 9{x) is the step 
function, f{e) = 1/(1 + exp(£:/T)) is Fermi function, T is temperature, Ri is the normal 
tunnel resistance of ith junction, and the energy gain W is given bylll 



Here Ci and C2 are the junction capacitances, Cs = Ci + C2 is the total island capacitance, 
and Qq is the total background charge which accounts for the initial background charge 
Qoo and the charge induced by the gate voltage Vg, Qq = Qoo + VgCg (for definiteness we 
consider the gate capacitance Cg as being added to Ci, although an arbitrary distribution 
of Cg between Ci and C2 is possible in calculations^!) . 

At low temperatures the quasiparticle current in SSS SET appears only above the voltage 
threshold 



(The last equation is the condition Wj^{n) = 2A(T).) 

At temperatures T comparable to the critical temperature Tc, the concentration of the 
thermally excited quasiparticles becomes considerable, and this modifies the shape of the 
/-F curve, in particular, creating several additional features at < Vt. As an example. Fig. 
HI shows the theoretical I-V curve of the symmetrical SET-transistor with A(T) = l.Se^/Cs, 
T = OAe'^/Cj] for several values of Qq (relatively large ratio A(T)/ (e^/Cs) is chosen to show 
more features). One can see two types of features: peaks (marked by arrows down) and 
steps (arrows up). 

The peaks positions constitute two seriesS 





Vt = mm{V^/\Vt>4A{T)), 



(6) 






(8) 



which correspond to the condition W^{n) = (obviously, the condition W~{n) = gives 
the same set of voltages). For such tunneling with zero energy gain the singularities of the 
density of states of two electrodes match (remind that we consider the same A(T) in all 
electrodes), that leads to the increase of the tunneling rate Tf{n) and explains the name of 
SM peaks. In BCS theory Tf{n) is formally infinite at Wj^{n) = (logarithmic divergence). 
Although the current through SET-transistor remains finite being governed by the stationary 
master equation (|^), the divergence of F would lead to very high and narrow center of the 
SM peak. To take into account the inevitable smoothing of the singularity of p{e), in Fig. 
m we assumed (phenomenologically) a small Gaussian inhomogeneous broadening of A(T) 
with dispersion w = O.OIA(T). The peak height depends very weakly (logarithmically) on 
w provided w ^ A(T). 

The origin of SM peaks is similar to that of well-known peaks! on the /- Fcurve of a single 
junction with different energy gaps Ai(T) and A2(r) of electrodes aiV = [Ai(T)— A2(r)]/e. 
In our case the energy gaps can be the same, and the energy shift is provided by the Coulomb 
blockade. However, this analogy is not complete. For example, in our case both singularities 
match simultaneously. Another difference is that the reverse process (tunneling back) also 
has a large rate, and the net transport is due to the tunneling through the other junction. 

The voltage position of SM peaks coincides with the position of the recently observed 
peaks0 in SSS SET at low temperatures when the parity-dependent current is due to the 
single quasiparticle created by the preceding tunneling event. 

At not too high temperatures the SM peaks are more pronounced within the voltage 
interval 2A(T)/e <V< 2A(T)/e + e/C-£ (see Fig. |l|). The lower bound is the condition that 
the tunneling through the other junction which restores the system into the initial charge 
state and gives the contribution to the net current, is sufficiently fast, W = eV > 2A(T). 
The upper bound is the condition that after this restoring the tunneling of the next electron 
through the same junction (which drives the system out of "resonance") has a small rate, 
W = eV — e^/Cj] < 2A(T). Hence, not more than two closely located peaks from the series 
given by Eq. (H) can be well pronounced on the /- V curve. 
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The important property of SM peaks is their specific temperature dependence. They 
should be absent at small T (because there are no thermally excited quasiparticles), and 
their height grows with T for some temperature range (see Fig. below) until they begin to 
decrease due to the suppression of superconductivity and/or correlation between tunneling 
events. Notice that the voltage position of SM peaks does not change with temperature 
despite the dependence A(T). 

One can see from Fig. [l| that the SM peaks are rather broad and have asymmetric shape 
so that they have longer and higher tail at the higher-voltage side. When the peak is not 
well-pronounced, this tail resembles plateau. When the SM feature is even weaker, it is seen 
small kink on the I-V curve (Fig. p. 

The other features seen in Fig. |l] are the step structures in the /- V curve which are 
similar to the step at V = Vt. Their positions satisfy the same condition Wj^{n) = 2A(T) 
and hence, the same Eq. (0) as for Vt. So the position of these two series of steps onV — Qo 
plane is just a continuation of the straight lines corresponding to Vt (they exist both above 
and below Vt). The steps in Fig. |I| are smoothed because of a finite w. 

Notice that while the steps corresponding to Eq. (0) are usually positive (increase of the 
current), they can also be negative - for example, when the step position is on the negative 
slope of an SM peak (the decrease of the current occurs because the charge state having the 
resonant tunneling rate becomes less probable). 

Besides the steps described by Eq. at relatively high temperatures the theory predicts 
an appearance of very small negative steps (both below and above Vt) at 

_(_!,(„ 4)1. (,) 

that corresponds to the condition W^{n) = — 2A(T). The steps are negative because the 
dependence T(W) given by Eq. (|) has the step down ai W = — 2A(T). The effect is very 
weak because of the factor exp(— 4A(T)/T) and also because of the existence of the opposite 
effect due to the simultaneous threshold condition Wj^{n — (—1)*) = 2A(T) (negative steps 
have not been observed in our experiment). 
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The aluminum-based single-electron transistors were fabricated using the standard two- 
angle evaporation technique. The details of the fabrication are given in Ref.§. Figure 0a 
shows the experimental dependence of the current on the gate voltage for different bias 
voltages at T = 650 mK. The largest feature seen is the Figure is the onset of the fast quasi- 
particle tunneling at V > Vt (Eq. (|^)). (Actually, we see peaks because of the small current 
scale of the Fig. for larger bias voltages they transform into plateaus with sufficiently 
sharp edges.) The well-pronounced peaks along the straight lines intersecting the abscissa 
axis at = — 14mV and Vg = 32mV are JQP peaks.i (They are due to Josephson coupling 
and, hence, are outside of the scope of the present paper. The position of JQP peaks is 
givenii'0 by Eq. (P) without the term 1/2 inside parentheses.) The step structures can be 
seen along the lines (see Eq. (|^)) which are the continuation of the main threshold Vt{Vg) 
(they start from the large features due to Vt in the upper part of Fig. |^a). And finally, the 
SM peaks are represented as rather broad features along the strait lines approximately in the 
middle (theoretically exactly in the middle) between JQP lines which intersect the abscissa 
axis roughly at ~ lOmV and Vg ~ — 40mV. Small SM peaks have been observed even 
above Vf. The SM features along the line with negative slope are more pronounced than 
along the line with positive slope possibly because of the difference in junction resistances. 
Similar measurements made at T = 50 mK do not show SM features as well as additional 
step structures while JQP peaks remain well-pronounced at \^ ^ 0.65 mV. 

Figure 0d shows the numerically determined positions of the maxima of the I — Vg curves 
from Fig. ^ on — plane. From the straight lines corresponding to JQP peaks (solid 
lines) we determine the junction capacitances Ci = 183 ± 4 aF, C2 = 117 ± 3 aF. The gate 
capacitance Cg = 3.5 aF which determines the gate voltage period of 46 mV is included 
into Ci because Vg has been measured from the outer side of Ci. Notice that the bias 
voltage corresponding to the intersection of two JQP lines directly gives the charging energy, 
e/Cs = 0.53 mV. The minimal Vt of 0.80 mV is used to determine the superconducting gap 
A(T) = 0.20 meV (minimal Vt corresponds to the edge of the almost vertical curved lines 
in the upper part of Fig. ||b). 
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Dashed lines in Fig. ^3 show the theoretical position VfJ^ of SM peaks calculated from 
Eq. (|^). We see that experimental peaks are located at somewhat higher bias voltages. This 
can be explained by several reasons. First, SM feature has a rather smooth shape, and, 
hence, the addition of any current component which increases with bias voltage leads to 
the apparent shift of the maximum to higher voltages (we also checked numerically that 
relatively large inhomogeneous broadening of A(T) leads to a similar shift). Second, the 
additional contribution to the position shift in Fig. |b can occur because the peaks are 
determined as the maximum current point over Vg, not over V. (The Vg change which 
decreases VfJ^^ also weakens Coulomb blockade in the same junction, hence, increasing the 
"background" current and leading to the apparent shift of the maximum position.) Finally, 
the third possible explanation of the shift (which we believe is most likely) is due to the 
difference between A(T) in the island and leads. Then each SM peak should split in two 
(there is some experimental evidence of such a splitting which is slightly seen in Fig. ^ and 
Fig. ^). Numerical simulations show that the peak corresponding to higher bias voltage is 
more pronounced (see Fig. ^) while the lower peak is possibly too small to be represented 
in Fig. The difference about 0.02 meV between the energy gaps would be sufficient to 
explain the experimental deviation. 

The experimental temperature dependence of SM peaks on / — curves is shown in 
Fig. Ha. Besides two SM peaks per period of Vg one can see two steps and two JQP peaks 
(the heights of JQP peaks are different because the Josephson coupling in one junction has 
been suppressed^). We see that SM peaks as well as steps grow with temperature. Solid 
lines in Fig. show the corresponding theoretical I — Qq curves calculated without fitting 
parameters (JQP peaks are not included in the model). The total resistance i?s = 605 kQ 
is obtained from the I — V curve and R1/R2 = C2/C1 is assumed. The gap A(0) = 0.207 
meV is used to get A(T) = 0.2 mV at T = 650 mK. The small broadening w = 0.03A(0) 
of the gap was used to eliminate the unphysical divergence of F at the peak center, while 
we did not attempt to fit experimental SM peak height by w (larger w decreases the height, 
though the dependence is quite weak for w ^ 0.05A(T)). The good correspondence between 
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Figs. 1^ and is an additional proof that the observed peaks are really SM peaks. The 
dashed curve (corresponding to the top solid curve) illustrates the peak splitting due to 
different A(T) in the island and leads (difference of 20 meV is used). One can see that 
this assumption not only explains the peak position shift and traces of such a splitting in 
experiment, but also improves the agreement for the peak height. 

Let us mention that the height of thermally activated JQP peak (at V < 2A/e + e/2Cs) 
as a function of V {Vg is varied correspondingly) should also exhibit the SM feature at 
V = e/2Cj] because at this voltage the tunneling of the second quasiparticle in JQP cycle 
is at resonance. There are some traces of such an increase in Fig. |^a and also one can see 
qualitatively that in Fig. ||a the JQP peaks start to decrease crudely atV< e/2CE (because 
of the thermal broadening of SM feature this smooth boundary moves to lower voltages by 
6V ~ TCY./Ci ~ 0.1 mV). Similar behavior should be expected for the height of thermally 
activated steps with SM feature at V = 2A(T). 

In conclusion, we observed SM peaks on I — Vg dependence of SSS SET at temperatures 
comparable to T^. The shape and position of the features agree well with the theoretical 
prediction. 

We thank the group from the University of Jyvaskyla for the information about their 
experimental resultslii prior to publication. The work was partially performed under the 
management of FED as a part of the MITI R&D Program Superconducting Electron Devices 
Project supported by NEDO. 
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FIGURES 

FIG. 1. Theoretical I-V curves of the symmetrical SSS SET with A(r) = l.Se^/Cs, 
T = 0.4e^/Cs atV<Vt for several Qq taking into account only quasiparticle tunneling. Notice the 
presence of SM peaks (marked by arrows down) and steps (arrows up). Small phenomenological 
smearing w = O.OlA(r) of the superconducting gap is assumed. Curves are shifted vertically for 
clarity. 

FIG. 2. (a) - The experimental dependence of the current / on the gate voltage Vg for SSS 
SET at r = 650 mK. The bias voltages V range from to 0.828 mV with the step of 7.08 /xV. The 
curves are shifted vertically by A7(pA)= 281 x V{inV). (b) - The positions of the current maxima 
on V — Vg plane. The solid lines fit JQP peaks. The dashed lines show theoretical position of SM 
peaks. 

FIG. 3. (a) - Experimental I — Vg curves for V = 0.69 mV at different temperatures which 
show two SM peaks (and also two JQP peaks and two steps) per period. The temperatures in 
mK (from top to bottom) are: 712, 684, 640, 605, 571, 532, 495, 462, 426, 386, 345, 303, 97. 
Notice that the height of SM peaks and steps grows with temperature, (b) - Solid lines show the 
corresponding theoretical prediction without fitting parameters. The JQP peaks have not been 
included in simulations. Small smearing w = 0.03A(0) is assumed. Dashed line illustrates the 
peak splitting due to different A(r) in the island and leads (±10 meV difference is used). 
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